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We examine quantum statistics of optical photons emitted from atomic ensembles which are 
classically driven and simultaneously coupled to a two-level atom via microwave photon exchange. 
Quantum statistics and correlations are analyzed by calculating second order coherence degree, von 
Neumann entropy, spin squeezing for multi-particle entanglement, as well as genuine two and three- 
mode entanglement parameters for steady state and non-equilibrium situations. Coherent transfer of 
population between the radiation modes and quantum state mapping between the two-level atom and 
the optical modes are discussed. A potential experimental realization of the theoretical results in a 
superconducting coplanar waveguide resonator containing diamond crystals with Nitrogen- Vacancy 
color centers and a superconducting artificial two- level atom is discussed. 
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I. INTRODUCTION 

Primary requirements of successful quantum informa- 
tion technology, in particular of large scale distributed 
quantum networking, is to have efficient and reliable 
means of storage [1], processing [2[ and dissemination 
[3| of quantum resources. These tasks can most likely 
be solved in different modules, designed and optimized 
for their ideal operation. The composite quantum infor- 
mation device would be then a hybrid structure of such 
sub-systems |4|. As the individual modules can have dif- 
ferent time and energy scales of operation, the fundamen- 
tal challenge is to merge them together congruously. 

In this paper we consider an atomic three-level ensem- 
ble, with the so called A- type level scheme p, [6j. We 
argue that such a system can be a promising set up for 
controllable exchange of quantum resources between en- 
ergetically remote modules in a hybrid quantum infor- 
mation device. Basic idea is to exploit the presence of 
an extra energetic component in the A-scheme to bridge 
the large energy gap between the other two components. 
In particular we envision a typical case where one com- 
ponent is a two-level atom, corresponding to the station- 
ary quantum hardware and quantum information unit 
(qubit). It is coupled to the A system by a cavity field in 
the microwave domain. The other component could be 
a remote quantum communication node such as a satel- 
lite or quantum repeater along an optical fiber so that 
these require coupling to the A-system in the optical 
domain. Hence by choosing the extra available transi- 
tion also in the optical domain, we can consolidate these 
two domains. The A system could play a double faceted 
role. It can serve as the quantum memory in terms of 
the pure or, together with the two level atom, hybrid 
spin ensemble qubits |7|. In addition, it becomes a medi- 
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ator of communication between the stationary quantum 
hardware and the quantum communication channels for 
optical flying qubits interconnecting remote nodes of a 
quantum network. 

Due to dipole selection rules, A-type transitions can 
happen only in certain systems, either lacking, or with 
explicitly broken, inversion symmetry, such as artificial 
atoms [5|, Q, semiconductors J9l4ll| or chiral molecules 
[12l [l3j . Alternatively, one of the transitions can be of 
magnetic dipole or electric quadrupole type in the A- 
scheme J14J-KL6J. Such weak transitions are collectively 
enhanced for strong coupling due to the ensemble p}. 
We specifically consider the latter case of magnetic dipole 
assisted A-scheme, though our results would be valid in 
other cases, too. 

Some quantum statistical properties of single isolated 
A-type atomic ensembles have been explored recently 
m 03 • Coupling of atomic ensembles, in particular, 
three-level A-type ensembles to a two-level atom has been 
discussed in the context of quantum state mapping from 
quantum memories to charge qubit [18| in transmission 
line resonators jl9J-[22j. In this paper we examine quan- 
tum coherence and correlations of the optical photons 
emitted by a A-type ensemble coupled to a two-level 
atom in a microwave cavity. We also investigate the 
case of two A-type ensembles. Equilibrium and non- 
equilibrium situations are separately examined. Con- 
ditions of photon antibunching, particle entanglement 
|23l l24j , genuine bipartite and tri-partite mode entan- 
glement, W state [25] are revealed. Furthermore, quan- 
tum state mapping between the two-level atom and op- 
tical modes, as well as coherent transfer of population 
between the radiation modes are found. 

This paper is organized as follows. In Sec. [II] we de- 
scribe the model system and by eliminating the collec- 
tive atomic states, we derive the effective Hamiltonians 
for two- level atom and radiation modes. In Sec. Hill we 
present our results and discussions for the single and two- 
ensemble cases in the steady state and non-equilibrium 
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FIG. 1. (Color online) Identical three- level atoms in the en- 
semble interact with two quantized (indicated by a, b) and one 
classical (indicated by Q) electromagnetic fields in A scheme. 
The microwave photons further couple the atomic ensemble to 
a single two-level atom. The classical drives on C-A transition 
and on the cavity field are not shown. 



situations. In Sec. IIV( a potential experimental system 
for physical implementation of our results is described. 
Specifically, we consider an extension of recently real- 
ized Nitrogen- Vacancy (N-V) centers in diamond crystals 
strongly coupled to superconducting coplanar waveguide 
resonators (CPWG) |26j. Similar electron spin ensemble 
to on-chip superconducting cavity couplings, specifically 
coupling of Cr 3+ spins in ruby and N substitution (PI) 
centers in diamond, are also observed in parallel efforts 
|27j . Finally in Sec. El we conclude. 



II. A-TYPE ATOMIC ENSEMBLE COUPLED 
TO A TWO-LEVEL ATOM 

A. Model Hamiltonian 

We consider a collection of N identical three-level 
atoms interacting with three electromagnetic fields in A- 
type transition scheme delineated in Fig. [TJ The atomic 
ensemble is further coupled to a two-level atom through 
microwave photon exchange. The lowest level | C) is cou- 
pled to the intermediate level | B) by a magnetic dipole 
transition of rate <#,, while the transitions from the high- 
est level | A) to the lower doublet | B) and | C) are optical 
electric-dipole transitions of rates Vt and g a , respectively. 
We further suppose that optical fields are traveling waves 
while the microwave field is a standing wave of a cavity 
that contains the two-level atom and the atomic ensem- 
ble. 

Our model can be conceived as a modest generaliza- 
tion of the case of a A-type atomic ensemble, resonantly 
coupled to the optical fields, which is studied for the bi- 
partite entanglement of the quantized fields [6| . On the 
other hand, such a seemingly innocent addition of two 
level atom generalizes the system from two interacting 
linear oscillators to the case of a linear oscillator cou- 



pled to a nonlinear oscillator, and as such can cause re- 
markable changes in the quantum statistical properties 
of the fields. In comparison to A-type ensemble coupled 
to two- level atom [28[, A-type scheme allows for effec- 
tive interaction between a, b modes. In addition, we take 
into account additional external drives and consider the 
system as open with all the decay channels due to cav- 
ity losses or spontaneous emission decays, which causes 
subtle effects in particular in the steady state. 

The system can be described by a model Hamiltonian 
H = Hq + Hi + H p where (in units of K) 

H 0=Y1 W afe4 a fc + Wfc6t& + Y^ + 51 "^ ' (1) 

k x 3 

Hi = gatf + n y^ R^ e ^- f '-""*'> 

3 

+ E 9 a ^R%a + g b E R%h + h.c. (2) 
H p = Q A Y^ iJ^e i(fc ^'^-° At) + E h e~ iGjbt b + h.c. (3) 



Here, x = A,Bj = 1,..,JV, uJ a ^uj b ,uj A ,ojB,oJo are re- 
spectively the frequencies of the electromagnetic fields 
and atomic transitions relative to the lowest level | C). 
Probe field frequencies are denoted by Co a and 0)5. Ef- 
fective drive strength of the microwave cavity and the 
Rabi frequency of the drive for the \C) to \A) transition 
are represented by E^ and Qa, respectively. Population 
and transition operators of single atoms in the ensem- 
ble are respectively represented by Rxx and R X y , x ^ y 
with x, y = A, B, C. Optical and microwave photon an- 
nihilation (creation) operators denoted by ag, b (at, fet), 
respectively. The coupling coefficients are assumed to 
be spatially homogeneous. We suppose that the wave 
vectors of the optical fields satisfy the mode matching 
condition [6] k — kd = &a so that quantized optical field 
can be treated as single mode. More detailed conditions 
on three-dimensional analysis of mode matching between 
free space light and the atomic ensemble will be ignored 
for simplicity [29| • Suppressing the constant global phase 
factor by [30] | A) — >> ex.p(ikd • rj) \ A), we drop the k 
index from our notations and further considerations. 

There would be additional Langevin noise terms in the 
Hamiltonian, corresponding to the coupling to the con- 
tinuum of other modes. These will be taken into account 
effectively in the master equation formalism and shall not 
be shown here. 



B. Quasi-spin wave picture of the model system 

It is illuminating to portray the model system in 
terms of collective excitation operators [3l| A = 
X C a/VN,B = X CB /\fN,T = X BA , where X xy = 
J2j=i Rxy, is analogous to the Hubbard operators [32| in 
strongly correlated systems. In the ensemble, the set of 



Hubbard operators represent SU(3) symmetry of a spin- 
1 system. They can be associated with the Gell-Mann 
fundamental representation of the SU(3) Lie algebra by 
choosing U = Xca-,V = Xbc and T = Xba, so that 
Z7, V,T form the three spin-1/2 subgroups. Each obeys 
the SU(2) algebra with [L,tf] = -2L 3 ,L = U,V,T. 
They are interdependent by the commutations such as 
[W, U^\ = —T^. Analogous to excitons and magnons, 
SU(2) algebras of U and V sub-spins can be contracted to 
Weyl-Heisenberg algebra h^ to associate collective excita- 
tions with quasi-spin wave quanta |3l| . The system space 
becomes a semidirect product Lie algebra, SU(2)®h2. In 
the limit of low excitations to | A) and | B), relative to TV, 
Vs ~ N/2 and Us « —N/2. With the contraction param- 
eter 1/^/N, the operators A = U/y/N and B = V ] /\fN 
become approximately bosonic. They become indepen- 
dent, as the isospin subgroup is weakly populated. The 
isospin is consistently represented by T = B^A. 

Denoting the phases of the g,g a ,#5, fi, QA,E b by 
Og,Oa,QbiQn,QA,0bi respectively, they can be suppressed 



by b 



b ex.p (iO b ) , a 



aexp[i(O g - b )],B 



— > 



£>exp [— i(0 b + 6b)] , a — »• aexp [i(6 a + 6a)] and A — >• 
Aexp(— %6a)- Employing a unitary transformation by 
U = exp(-iu d t) with H -► WHU - iU^d t U, the Hamil- 
tonian becomes 



H = J2 ™xX ] x + uj b b ] b + uj B B ] B + ^-cr z , (4) 

x=a,A 

H x = gatf + g a A ] a + g h B*b + fte^A^B + /i.e., (5) 

H p = Q A e-^ At A^ + E h e- iC ° bt b + /i.e. (6) 



where the coupling coefficients g a ,b ^ \Ng a $ are collec- 
tively enhanced; (j) = Qn + A -O b - Q b , m a , A = u a ,A - oJd 
and tz7^ = (I>a — uud. Even though it can play curious 
role on quantum correlations [6], for simplicity we take 
4> = 0. We further assume three-photon resonance con- 
dition voa — &a —Ud = oJb f° r the classical driving fields. 



C. Effective Hamiltonian for the model system 

A-type ensemble can induce an effective photon - pho- 
ton interaction mediated by the quasi-spin wave back- 
ground [6|. In order to show that explicitly, let us intro- 
duce the normal modes p± = u±A + v±B by the Bogoli- 
ubov transformation [33] with the relations U- = — v+ 
and V- = ix + , where u± = ^/{rj ± (wa — ujb))/2t] and 

T]= > /(tx7 A -£J B ) 2 +4n 2 . 

The model Hamiltonian can be rewritten as H = Hq + 



H x 



H p where 



H = w a a ] a + cj^b + ^ ^aP^a + y^, (7) 



A=± 



#i = gafr 1 " + ^2 P\(9a\a + #6a6) + h.c, (8) 



A=± 



H P =Y^ n A xe~^ At px + Ebe-^b + /i.e. (9) 

A=± 



Here ft± = (zu A +^b ±^)/2, #aA = 9 a ux,9b\ = 9bV\ and 

^AA = ^4^A- 

Effective coupling between the electromagnetic modes 
can be obtained by a unitary transformation to elimi- 
nate the interactions between the electromagnetic and 
quasi-spin- wave modes. We introduce an anti-Hermitian 
operator S 



S =Y1 ( Xxa + y xh )A ~ h ' c 'i 



(10) 



A=± 



where the unknown coefficients x\ and y\ are to be de- 
termined by the Frochlich - Nakajima canonical transfor- 
mation condition 133-1361 



[Hq,S] = -^2 p{(ga\a + gb\b) 



h.c, 



(ii) 



A=± 



which yields x \ = g a \/(^a-^x) and y x = g b \/(uj b -Q\). 
Assuming ft ^> g a: g b: an effective Hamiltonian H e ff = 
exp (-S)Hexp (S) can be written up to the second order 
in x x , yx, 

H e ff = w' a a)a + u' h b^b + — -a z + g(crtf + h.c.) 
+ (E' h e- i(:jbt b + £; a e-^ A ^a t - Ja f 6 + h.c) 



A#a* 



^(GxcTpi + Epxe-^pl + h.c.) 



(12) 



with A, /i = ±. The coefficients in i? e // are listed in 
the appendix [A[ The indirect coupling between the 
quasi-spin-wave modes and the electromagnetic modes 
is through their common interaction with the two-level 
atom. In the large Q limit, the quasi-spin waves are 
strongly off-resonant with the two-level atom and practi- 
cally uncoupled. More formally, in the dispersive regime, 
Q — ujo ~ Q ^> g, employing another Frochlich trans- 
formation, such an interaction can be reduced to a Stark 
shift of two- level atom resonance proportional to the pop- 
ulation of the quasi-spin normal modes, that can be ne- 
glected for weak excitations. Influence of such sequential 
canonical transformation on radiation modes would be 
negligible provided that Q ^> g a ^> g bl g. In addition, 
effective weak drive on the two-level atom would be neg- 
ligible. The effective Hamiltonian, keeping the leading 
terms only, is reduced to 



H e ff = (m a - £j b )a^a + (u b - cu b )tfb + 

+#(crfr t + h.c.) - J(a ] b + h.c.) 
+E a (a^ + a)+E b (tf + b). 



ujo -Wb 



(13) 



Here, three-photon resonance condition for the driving 
classical fields is employed and the Hamiltonian is written 
in the frame rotating at uo b . We further suppose that 
three-photon resonance holds for the radiation fields such 
that uj a = ujd + oo bl then w a — uo b = uj b — Cj b = S is 
defined. We find ujq — oj b = S — 5 b . According to Fig. [TJ 



S b = cjfo-cjo, A = oj d -uj AB , A a = uj a -uj A ,A b =uj h -u B , 
in terms of which we express the three-photon resonance 
condition A a — A — A5 = 0. 

The Hamiltonian of Eq. [13] has a zero-eigenvalue so- 
lution in the single excitation manifold, or the so-called 
dark state, which is given by |28j | 



DS) = (cos aa — sin aay \ vac), 



(14) 



with a = — arctan(g/J) Presence of such a dynami- 
cally fixed state has been exploited to propose an adi- 
abatic transfer protocol to map the quantum informa- 
tion stored in a molecular A type ensemble to a charge 
qubit of two-level atom, and vice versa. Our treatment 
of cyclic ensemble reveals a complement to this proto- 
col by allowing quantum information transfer from two- 
level atom stationary qubit to optical flying qubit. A 
state (x \g) -\-y |e)) |00) a5 can be adiabatically transferred 
into {x |0) a + y |l) a ) \g) |0) b . The decoherence channels 
due to two-level atom and the microwave cavity are sup- 
pressed. The effective loss channel of the optical photon, 
transferred from the spontaneous emission of the level 
\A) during the Frochlich transform, is compensated with 
the effective drive E a . Taking into account time depen- 
dent interaction coefficients and time-dependent canoni- 
cal transformations, an optimal fidelity of the state trans- 
fer can be characterized. 

Alternative to the derivation based upon sequential 
Frochlich transformations as presented here, it is pos- 
sible to get the effective Hamiltonian with a single step 
Frochlich transformation or by reducing a general four- 
color (a, 6, A, B) polariton to a two color (a, b) one [6[ for 
strong drive Q conditions. It is straightforward to gener- 
alize this model by either of these methods to the case of 
a pair of cyclic atomic ensembles in the microwave cavity, 
coupled to the same two-level atom. We find 

H *ff = J2 W aj a ] a 3+u' b tfb+^-(T z 
i=l,2 ' Z 

+g(ab t + h.c.) - ^ Jj{a]b + h.c.) 

j=l,2 

+ J2 E M + a j) + E ^ + h )- ( 15 ) 

i=i,2 

For identical symmetric ensembles, w a j = vo ai Ej = E a 
and Jj = J. 

The single ensemble effective model reduces to the 
Jaynes-Cummings (JC) model for the optical mode in 
the dispersive limit of the microwave cavity. Two- 
ensemble model becomes a two-site Bose-Hubbard (BH) 
model, with pure and cross-Kerr interactions of the type 
(a\ai) 2 ,i = 1,2 and ajaia^Gte- A common proposal to 
realize a BH dimer is to use JC-BH systems of tunnel 
coupled nonlinear cavities, where the nonlinear ity is pro- 
vided by a two- level atom in each of them J37|, [38j . In 
contrast to such direct methods to realize atom-optical 
analogs of Josephson networks, here there is only one 
two-level atom providing shared nonlinearity by each op- 



tical mode, which simplifies architecture of identical ho- 
mogeneous coupled subsystems. 

BH dimer can be described as a pseudo-spin sys- 
tem with quadratic nonlinearities using pseudo-spin 
Schwinger operators S+ = a[a2,S z = (a\a\ — a\a2)j2. 
When sufficiently strong nonlinearity present, it leads to 
multi-particle entanglement between the distant optical 
modes via the so-called axis twisting spin squeezing route 
[23l . l24{ . The question of whether multi-particle entan- 
glement can still be found in the non-dispersive regime 
is not as obvious to answer due to lack of immediate 
mapping onto a spin model. We shall explore that case 
numerically in the following section, together with other 
types of quantum correlations as well as with quantum 
statistical properties of radiation modes in the effective 
models. 

Generalization of these models beyond two-ensembles 
is in principle possible. However, the quasi-spin wave 
mediated effective tunnel coupling here is long range be- 
yond nearest neighbors and hence direct relation to sim- 
ple Bose-Hubbard models is not immediate. 



III. RESULTS AND DISCUSSIONS 

Assuming the effective model remains Markovian and 
Born-Markov approximation is applicable, dynamics can 
be investigated by a quantum master equation 



i[H,p] + J2^ V W 



x\p + <yV[*]p+?±V[a z ]p(W) 



Here n x is the photon loss rate of the mode x = a,b, 
71 and 7^ are respectively the relaxation and pure de- 
phasing rates of the two-level atom. Dissipation terms 
V[x]p = (2xpx* — x^xp — px^x)/2 are the Liouvillian su- 
peroperators for damping in Lindblad form. 

Frochlich transformations used in the derivations of the 
effective models influence the decay rates. In the lead- 
ing order, we estimate K a ~ 7^(^ a /^) 2 , where ja is the 
excited state decay rate of the A-type atom. We choose 
g a /£L = 0.1 so that n a = O.OIja- We assume n a = n^. 
In addition, we suppose E a = E b = 0.1^5, which corre- 
sponds to weak driving condition for which empty cavity 
photon number is 0.01 on resonance. Such strict equali- 
ties are to reduce the number of independent parameters 
in the simulations but not crucial in practice. It is suf- 
ficient to arrange field-atom detunings and the driving 
field strengths (The formulae in the Appendix [A[ can be 
referred) to ensure sufficient effective drive E a can be 
translated from A to a so that the optical radiation can 
survive in the steady state against the decay channels. 

We numerically solve the master equation using the 
"Quantum Optics Toolbox" [39] and determine the 
steady state and dynamical behavior of the system. Our 
typical results will be reported below for the cases of sin- 
gle and two-atomic ensembles separately. For all cases 
we fix J/2tt = 1 MHz. 



A. Single Atomic Ensemble 

1. Steady State Coherence and Correlations 

Taking S h = J, K h = K a = (2tt)0.4 MHz, 7i/2tt = 0.02 
MHz, 70/27T = 0.3 MHz, steady state solution of the 
master equation is determined, in a truncated Fock space 
of 3 photons in each mode, for a range of g/2ir G {0, 4} 
MHz and S/2tt G {-5, 5} MHz. 

Populations of the radiation modes are shown in Fig. [2] 
in the (£, g) phase space. Dashed lines indicate the probe 
resonances with the single excitation manifold. The rel- 
atively low transmission along a particular resonance is 
due to our choice of equal drive strengths. For a simple 
illustration, let us consider S^ = for which 
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are the dressed atom-polariton operators [40j of the single 
excitation manifold. For our choice of S b = 0, q D \vac) is 
the dark state of the manifold. The states generated by 
the q± are with ±K energies where K = ^/ g 2 + J 2 . For 
g < J we find a = (g + + g_)/v // 2 and b = (q+ + q-)/y/2 so 
that the combined drive terms for the a, b modes becomes 
y/2E a {qj r + q+) if E a = E^. In that case neither q- nor 
qo are driven and hence they cannot survive to steady 
state. Fig. [2] confirms with this argument as it exhibits 
that only bright polar iton branch is the upper one, while 
the lower two remains in the dark. With the increase 
of #, a mode drives mainly the qo polariton. This is 
seen in the Fig. [2] as the bright middle polariton branch. 
The upper branch of q+ is more strongly driven relative 
to q- and correspondingly associated with the bright b 
line in the middle figure of Fig. [2j One can selectively 
populate different polariton branches, for example taking 
Eb = —E a yields bright b mode along the lower polariton 
branch. 

The right figure in Fig. [2] is the behavior of popu- 
lation fraction defined to be w = (n a — nb)/(n a + n&), 
with n x = (x^x),x = a,b. As long as g < J there is 
no distinct radiation phases in the (5, g) space. When 
g > J however, radiation in the system is dominated by 
the microwave mode along the upper and lower polariton 
branches. In between, there is a wedge like region about 
the middle polariton branch, where the radiation is dom- 
inated by the optical mode. The main effect of non-zero 
St, would be to change the slopes of the probe resonance 
lines, and as such the widths of the optical and microwave 
phases, but the above descriptions would be valid. 

We characterize quantum statistical coherence proper- 
ties of the radiation modes in the steady state in terms of 
the second order coherence function at zero time delay, 
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2. (Color online) Steady state populations n a (Left Fig- 
rib (Middle Figure), and their relative fractional pop- 
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FIG. 

ure) , 

ulation n a b = (n a — nb)/(n a + rib) (Right Figure) of the 
optical and microwave modes at J/2tv = 1 MHz, Sb = J, 
Kb = Ka / = (2tt)0.4 MHz, 7i/2tt = 0.02 MHz, 7^/2^ = 0.3 
MHz. Fock spaces for the radiation modes are taken to be 
four dimensional. Dashed lines indicate the probe resonances 
with the single excitation manifold. g/27T and S/2tt are in 
MHz units. 



the von Neumann entanglement entropy [2j, 

S = -Tr (p a log 2 p a ), (20) 

and the genuine two- mode entanglement parameter [41] , 

Kb =| (a^b) | 2 -<aW&>. (21) 

Here, p a = Tri ) (Tr (T (p ss )) is the reduced density matrix 
of optical mode a evaluated by tracing out atomic (a) 
and microwave (b) degrees of freedom. S is an entangle- 
ment measure for pure bipartite states. As we have three 
sub-systems, we can interpret it as a measure of impurity 
in the steady state, analogous to linear entropy. On the 
other hand, we find that the two-level atom is weakly 
excited everywhere in the J, g phase space for the steady 
state. Thus, S bears the signs of profound quantum cor- 
relations between the radiation modes, especially in the 
g < J regime where it should be a reliable measure of 
mode entanglement. A more clear signature of genuine 
two- mode entanglement can be found when A a 6 > 0. 

Mode entanglement is due to correlations in the occu- 
pation number space in the second quantization frame- 
work. In the first quantization or in the particle picture, 
scattering of particles can lead to multi-particle entan- 
glement. Though we do not have such direct particle 
interactions for photons, effective interactions that lead 
to photon blockade effect [42|, |43| can yield similar multi- 
particle entanglement of photons. We want to explore if 
increased particle correlations among the photons can be 
found with the increasing anharmonicity parameter g in 
regimes of sub-Poisson statistics. To characterize particle 



entanglement, we use spin squeezing parameter [23l l24j , 
generalized to the case of dissipation, [44] defined to be 
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(22) 



When ^ 2 < 1 the state is said to be particle entangled, 
so that it has non-separable density matrix in the first 
quantization. Here, the pseudo-spin operators J x ,y,z are 
the usual Schwinger representation for the two bosonic 
modes a and 6, 



J x = -(a+fc + a&t), 

J z = l -((Ja-tfb). 



(23) 
(24) 
(25) 



(N) = (a^a + tfb) is the total number operator subject 
to losses. The squeezing is in general defined for an arbi- 
trary axis, which makes it possible to optimize. However 
it is more practical to examine a particular axis in the 
lab frame, which we choose it to be the x-axis. 

According to Figj3j S is strongest about Mq = (g,5) ~ 
(J/2, J). For g > J, it decreases, but occur predomi- 
nantly along the polar iton branches. The decrease with 
g is due to reduction of the Hilbert space dimensionality. 
All the three subsystems, two- level atom and the radi- 
ation modes, are excited About Mo point which can be 
seen in Figs. [2](3j Three component atom-cavity "molec- 
ular" structure reduces to two-component forms as g in- 
crease. The middle and right polariton branches becomes 
dominated by a and b modes, respectively. Two- level 
atom excitation is mainly along the right branch. The 
middle figure indicates genuine two-mode correlations are 
weakest about Mq and stronger along the middle branch,. 
The opposite behavior of A a 6 and S confirms the strong 
influence of atomic component on S. X a b < reflects that 
a — b correlations are atom mediated rather than being 
genuine. Though not visible in the scale of the figure, 
A a & > is found in the vicinity of (S = l,g = 0) point, 
for which decay channels of the two-level atom is turned 
off and true two-mode entanglement is survived, though 
barely, in the steady state. 

The second order coherence functions for the modes 
a, b together with the spin squeezing parameter of multi- 
partite entanglement are shown in Fig. 2J We plot 
-log(^(0)),-log(^(0)), and -log(g) so that their 
positive values respectively mean sub-Poisson statistics 
for modes a, b and particle entanglement of optical and 
microwave photons. 

The behavior of the g 2 (0) for g > J is typical in JC 
model. For g < J, g 2 (0) - gl(0) < 0. The equiva- 
lence of the g 2 (t) to Mandel Q parameter [45| associates 
anti-bunched photons with g 2 < with sub-Poissonian 
statistics. The only significantly populated anti-bunching 
region for mode a is about S ~ —0.5, at the beginning of 
the optical domain. The other two anti-bunching regions 
of the optical mode are along the microwave domains and 
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FIG. 3. (Color online) Same as Fig. [21 but for the von Neu- 
mann entanglement entropy S (Left Figure), genuine two- 
mode entanglement parameter X a b (Middle Figure), and the 
excited state population a ee — cr^a (Right Figure). X a b is 
multiplied by 10 5 and a ee is multiplied by 10 3 for visibility. 
Dashed lines indicate the probe resonances with the single 
excitation manifold. g/27r and S/2tt are in MHz units. 
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FIG. 4. (Color online) Same as Fig. [21 but for second-order 
coherence functions g^ifi) (Left Figure), ^(0) (Middle Fig- 
ure), and the spin squeezing parameter £ x (Right Figure). 
These are plotted as — log(pa(0)), — log(^(0)), and — log(^). 
g/27i and 5/2tt are in MHz units. 



hence relatively weakly populated. The coherence prop- 
erties of the microwave photon is transferred to the opti- 
cal one via their common interaction with the two-level. 
Collective enhancement of the interaction coefficients in- 
duces bosonic, linear oscillator, character to the A-type 
ensemble, which is harmful for sub-Poisson photon statis- 
tics. By addition of the two-level atom, dressed system 
regains its lost anharmonic character. Enhancement of 



both J and g optimizes the anti-bunching of the opti- 
cal photons. Despite relatively small coupling constants, 
the minimum ^(0) ~ 0.8 and #&(0) ~ 0.1 are within 
the same range with two- level nonlinear oscillators. A 
promising route to enhance these numbers could be to 
consider a more efficient source of anharmonicity, such 
as the TV- type 4- level system [46]. 

Nonlinearity of the JCM is associated with the anhar- 
monic spacing of the energy levels. One prominent effect 
of anharmonic level spacing is photon blockade [42| which 
is the suppression of the sequential absorption of two pho- 
tons into the system. In our case the anharmonicity influ- 
ences the normal modes of the coupled a, b modes similar 
to a case of photon blockade in a two-mode cavity. The 
degree of second order coherence is equivalent to M andl 
Q parameter which is a relative measure of number fluc- 
tuations. Photon blockade is associated with sub-Poisson 
statistics of the anti-bunched radiation. Since the normal 
mode transformation is a rotation in the mode space, sub- 
Poisson statistics that is found in a normal mode can be 
translated into the original modes. Such quantum noise 
description of the effect allows us to appreciate the local- 
ization of the relatively strong particle correlations about 
the polariton branches. 

Due to strong dissipation channels and relatively weak 
coupling strengths the steady state exhibits low level of 
quantum correlations among the subsystems. On the 
other hand, some of its promising features can be found 
in early time quantum dynamics of the system. Our ob- 
jective in the next section would be to go beyond steady 
state and explore some non-equilibrium cases for our 
model. 



2. Non-Equilibrium Coherence and Correlations 

We investigate dynamics of quantum correlations 
among the radiation modes in the single ensemble case for 
initial state |^(0)) = |11#). Truncating the Fock space 
for each radiation mode at 7 photons, quantum master 
equation is solved by the direct integration. 

To comprehend clearly the effect of g on the correlation 
dynamics, we first consider the case /%, ^5, 7, 7^ = 0, and 
accordingly, E a ,b — 0. We choose £5 = for which the 
probe resonance with the upper polariton branch is at 
5-. 
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By calculating the <^, we find that there is no multi- 
particle entanglement for g < J. We examine ^ for g > 
J in Fig. [5] We plot — log 2 £% so that any positive value 
of it would be the indicator of the particle entanglement. 
Strongest multi-particle correlations exist when g ~ J. 
Further increase of g gives weaker squeezing. While the 
entanglement happens earlier at larger g, its duration 
severely decreases with g. 

In the regime g < J lacking multi-particle entangle- 
ment, there is strong mode entanglement as depicted 
in the first column of Fig. [6j where we plot the S 
together with the fidelities Fj = Tr(\ipj) (ipj\ p) of the 




FIG. 5. (Color online) Time and g dependence of ^ when 
there is no dissipation in the syste m for a n initial state 
|^(0)) = |ll) ab \g) at 6 b = and S = ^Jg 2 + J 2 with J/2tt = 1 
MHz. In the figure g/2n is in MHz and t is in jis units. 



states \i/Jj=o.a) = \Ug) , |10e) , \20g) , \02g) , |01e), respec- 
tively. For g ~ there are two maximally entangled 
two-particle states [6|. The one appearing at the peaks 
S ~ 1.5 is \i(j) = (\20g) + \02g) ± i \llg))/y/S. It is en- 
tangled in three dimensional space so that the maximum 
S = /o#2(3) ~ 1.5. While the one at the peaks S ~ 1 is 
of the form \ip) = (\20g) + \02g))/y/2 and is entangled in 
two dimensional state space. At regular times the sys- 
tem is disentangled to |^o)- Strong and long-lived mode 
entanglement can be found when g ~ J as shown in sec- 
ond column of the Fig. [6] Splitting of the degenerate 
states helps to preserve three dimensional structure of 
the subspace of entanglement to ensure highest available 
S ~ 1.5. In addition, destruction of the disentanglement 
ensures high S can be maintained relatively longer dura- 
tion in comparison to g ~ regime. For g ^> J, accord- 
ing to the last column of Fig. [6j dynamical behavior of 
S resembles qualitatively g ~ regime, but for the lower 
S ~ 1 values. The states at these points are constructed 
by I llg) and |01e) for the middle, lower peak, and by 
I lOe) and \02g) for the higher peaks. They are both in 
two dimensional state spaces so that S ~ 1. 

We now examine the behavior of correlations in the 
presence of dissipation. A typical result for our simula- 
tions is shown in Fig. [7] indicating the dynamical behav- 
ior of S for the same parameters with the lossless case but 
for ft a /27r, ft6/27r,70/27r = 0.1 MHz. If we use the larger 
decay rates, same with those used in the steady state 
analysis, and take the same 5^ = J with the steady state 
case, we find the result in Fig. [8j We see that the fea- 
tures of mode entanglement discussed in dissipationless 
case is vastly available in the case of moderate damp- 
ing. In the larger damping case, except the early times, 
entanglement cannot be obtained. Some signatures of 
the prolonged entanglement duration is still found in the 
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FIG. 6. (Color online) Time dependence of S 

and fidelities Fj,j = 0,1,2,3,4 of the states 
\llg), |10e) , 1 20^) , 1 02^) , |01e), respectively. There is 

no dissipation in the system an d the i nitial state is 
|^(0)) = |ll) ob \g). 5 b = and S = y/g 2 + J 2 with J/2tt = 1 
MHz are used. In the figure t is in /xs units. The left, middle 
and right columns are for g/2-K — 0, 1, 10 MHz. Dash-Dotted 
curves in the first, second and third rows are for S, F2, F4. 




FIG. 8. (Color online) Same as Fig. [3 but for ft a /27r 
0.4, ftb/27r = 0.4, 70/2tt = 0.3 MHz and 5 b = J. 
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FIG. 9. (Color online) Same as Fig. [71 but for the squeezing 
parameter £j!. 



FIG. 7. (Color online) Time and g dependence of von 
Neumann entropy for an initial state 1^(0)) = |ll) a & \g) at 
S = y/g 2 -j- J 2 and £b = 0. In the figure g/27r is in MHz and t 
is in /xs units. The decay rates are ft a /27r, Kb/^TY, 7</>/27r = 0.1 
MHz. 



g ~ J even under such strong damping condition. 

Figs. [9l fim show that multi-particle entanglement can 
still be found at certain g values at particular time 
domains for weaker damping situation, while no spin 
squeezing and entanglement remain in the larger dissi- 
pation case. 



B. Two Atomic Ensembles 

We shall now consider the case of two identical, sym- 
metrically placed cyclic ensembles. In contrast to a — b 
coupling, there is no direct a\ — ci2 interaction now. We 
did not find multi-partite entanglement between the pho- 
tons of ai and (12 modes. We focus on mode entanglement 
and population dynamics in this case. Since the system 
has three mode radiation modes, we take into account 
possibility of genuine three- mode entanglement. We use 
an extended set of mode correlation parameters [41| 



A 16 = I (a\b) | 2 -(ajaifc+fc), 
Mb = I (a\b) | 2 -(4a 2 &t&), 
A12 = I (a|a 2 ) | 2 -(ajai4a 2 ). 



(26) 
(27) 
(28) 
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FIG. 10. (Color online) Same as Fig. [8] but for the squeezing 
parameter £ x 



When any two of them become simultaneously positive, 
genuine three-mode entanglement is found. Typical ex- 
ample is the so called three parity W state [25j which is 
defined to be W = (|100) + |010) + |001>)/a/3. In our 
case the Hilbert space is enlarged by the two-level atom. 
We denote basis states with |ni, 715,712 ?s) where rtx^^b 
are the occupation numbers of the radiation modes a\^b 
and s = e,g. 

We solve the density matrix equations using quantum 
trajectory method. We take 25 trajectories and trun- 
cate Fock space at 2 photons in each radiation mode 
ai } 2,6. We fix £5 = and S = \/ J 2 + g 2 for a given 
g for which the probes are in resonant with the upper 
polariton branch of the single excitation manifold. For 
arbitrary A there are four such branches now due to ad- 
ditional optical degree of freedom in comparison to the 
single-ensemble case with three polariton branches. We 
explore dynamical evolution of two initial preparations of 
the system: (i) |^(0)) = \010g) and (ii) |^(0)) = \100g). 

In the first case, our simulations reveal that when 
g < J population of the microwave cavity photon is sym- 
metrically split into to the optical modes as shown in Fig. 
ITT1 The effect is an optical analog of coherent population 
trapping [47} . At large g > J splitting is suppressed and 
the cavity photon is predominantly making localized in- 
teractions with the two-level atom. To comprehend this 
effect let us write the dark state \DS) in the single exci- 
tation manifold for the two-ensemble model as follows 



\DS) 



1 



VJ 2 



( j t 9 f 9 f 



\vaclW) 



with \vac) = 1 000g) being the vacuum state of the com- 
posite system. For g ^> J, \DS) becomes \DS) ~ 
-(\l00g) + \001g))/y/2, while for g <C J it gets \DS) « 
|000e). This suggests that in the regimes where the long 
range interaction J dominates over the short range in- 
teraction g, the dark state |000e) is inaccessible by the 



FIG. 11. (Color online) Time dependence of populations of 
the radiation modes m = (a\ai) (Left Figure), rib = (b^b) 
(Middle Figure), ri2 — (0^2) (Right Figure) for an initial 
state I ^(0) = |010p). We take Ji, 2 /2tt = 1 MHz, 5 b = 
0,(5 = y/g 2 + J 2 , «i, 2 ,6 = (2tt)0.4 MHz, 7 /2tt = 0.02 MHz, 
7^/2tt = 0.3 MHz, Ex = 0.1k x ,x= 1,2,6. 



initial state |^(0)) = \010g). As such, |^(0)) only yields 
a coherent transfer of b population to the a\^ modes. 
According to Fig. [TT] population transfer is complete in 
about t ~ 1 /is. Population exchange is a coherent pro- 
cess happening at regular time intervals of every 2 J. At 
lower damping rates, such as 7, ^1,2,5/^ ~ 0.1 MHz, we 
find a few more cycles of transfer. 

When the tunneling from b to a\^ is complete, the 
system radiates genuine two-mode entangled optical pho- 
tons which can be verified by the Fig. [T2I We veri- 
fied S ~ 1 is a reliable bi-partite entanglement signature 
at the corresponding times for low (7,^1,2,6/2^ ~ 0.1 
MHz) damping. During the course of transfer, at about 
t rsj 0.5 /is all three modes are populated. At such 
times, as the excited state of the two- level atom remains 
in the dark, three parity W state with genuine three- 
mode entanglement is found. Corresponding S ~ 0.7 
is relatively lower than the maximally entangled two- 
mode state. Due to its robustness against single parti- 
cle losses relative to Greenberger-Horne-Zeilinger (GHZ) 
states IH SI (\GHZ) = (|000) + |lll»/ v / 2) and its 
potential role in quantum networks, W states have ap- 
pealing features in quantum information science. They 
are proposed for bi-partite entanglement distillation [50] . 
Our results provide a promising setting to realize these 
goals. 

In the g > J case, \DS) « — (|100^> + \001g))/V2 
remains in the dark. The initial state is predominantly 
coupled to |000e). The radiation energy is trapped in 
the microwave cavity (b) mode and exchanged locally 
with the two-level atom. Associated Rabi oscillations are 
shown in the middle figure of Fig. [TTJ Mode entangle- 
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FIG. 12. (Color online) Same as Fig. [TT]but for the genuine 
mode entanglement parameters A12 (Left Figure), \2b (Middle 
Figure), and Ai& (Right Figure). 



FIG. 13. (Color online) Same as Fig. 111! but for an initial 
state I ^(0) =| 110) \g). 



merit parameters practically vanish in this case, as shown 
in Fig. [H 

We consider an asymmetric populated initial state to 
show that instead of two-mode entanglement of optical 
modes, optical and cavity mode entanglement can be em- 
phasized. For this aim we start with ^(O)) = |110g). 
Population dynamics at various g is depicted in Fig. [13l 

This localization of the cavity photon in the g > J Rabi 
oscillations regime can be used to block an optical mode 
to access the other distant optical mode and as such, the 
corresponding distant ensemble. The localizations of the 
radiation energies in g ^> J regime can be seen as block- 
ing of an optical mode interacting with an ensemble to 
access a distant ensemble and the corresponding optical 
mode via the microwave nonlinear cavity. In g <C J limit, 
the b mode still splits into two but its tunneling into a\ 
is balanced with the tunneling of a\ into b. There is no 
net population current between the a\ and b modes, so 
that we see no change of population of a\ mode in the 
left figure in Fig. [13j except its usual decay. Accordingly, 
A12 < and Aib < as indicated by the white regions 
about t ~ 0.5 /is in the left and the right figures of Fig. 
[141 While the middle figure demonstrates that there is 
genuine two-mode entanglement between the <22 and b 
modes at this time. In contrast to previous W state aris- 
ing at this time with three- mode entanglement, initial 
imbalance of the populations in the present case effec- 
tively removes one particle from the W state to distill 
or concentrate a two- mode entanglement out of it. At 
t — 1 /is the situation is similar to previous symmetric 
initial state and we see two-mode entanglement between 
the optical a\ and CL2 modes. The effect can be viewed as 
the presence of dynamically distinct and transferable bi- 
partite entanglement phases in the multi-mode system. 
While we consider here only the simplest symmetric and 
identical two-ensembles time dependent controllable dif- 
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FIG. 14. (Color online) Three-mode correlation parameters 
that can be used to reveal Genuine three-mode entanglement 
in the evolution of an initial state | ^(0) =| 110) | g). 



ferent interaction coefficients J\ , J 2 , g allows for more rich 
coherent information and population transfer protocols. 



IV. PROPOSAL OF EXPERIMENTAL 
IMPLEMENTATION 

In order to realize the model Hamiltonians discussed in 
the preceding sections, now we shall consider a straight- 
forward extension of a recent experiment demonstrating 
strong coupling of an N-V center ensemble to a coplanar 
waveguide resonator [26j. The setting is shown in Fig. [15j 
where two diamond crystals and a transmon qubit I5ll— 
[53j strongly coupled to microwave stripline cavity [19| . 
Diamond crystals containing color center defects serve 
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FIG. 15. (Color online) Transmon qubit and two diamond 
crystals couple to microwave stripline cavity of lenght L and 
width d. 



as the quantum memory, transmon qubit serves as the 
quantum hardware for rapid quantum information pro- 
cessing, and the microwave cavity photons serve as the 
local data bus. By interfacing this device with external 
optical fields, such stationary quantum information unit 
would have access to optical flying qubits for quantum 
networking and communication at free space communi- 
cation wavelengths. 

Strong coupling of a single diamond with Nitrogen- 
Vacancy (N-V) centers to a superconducting coplanar 
waveguide (CPWG) resonator has been shown recently 
[26j . We propose a a straight forward extension of the 
physical system of Ref. [26[ to the case of two diamond 
crystals. In order to place the diamond crystals to lo- 
cations where the magnetic field is strong and homo- 
geneous, and the transmon qubit to a similar location 
but for the electric field, a L = A cavity is considered 
suitable for our purpose. A ~ 50 mm is the microwave 
wavelength for the fundamental frequency of the cavity 
uji/2tt ~ 6 GHz. The crystals are assumed at the same 
size of 3 x 3 x 0.5 mm 3 . Canonically quantized magnetic 
field of the fundamental quasi-TEM mode of the cavity 
is given in the lab (cavity) frame by 



B(x, y, z) = B_\_(x, y)(b + 6' ) cos 



L 



(30) 



where B_\_(x, y) is the transverse field distribution. Place- 
ment of the crystals A/2 apart, at the magnetic field 
antinodes of the fundamental mode along the cavity axis 
(x), allows for sufficiently strong magnetic coupling to 
the cavity mode as well as adequate space for optical ac- 
cess to each crystal. At the same time, transmon qubit 
placed in the middle of the crystals would be at the elec- 
tric field antinode of the fundamental cavity mode. It 
is now well-established that, under certain conditions on 
the Cooper-pair box and circuit QED system, Jaynes- 



FIG. 16. (Color online) Energy level diagram of negatively 
charged N-V center. ZPL is the zero phonon line at 638 nm 
(1.945 eV) and ZFS is the zero field splitting at 2.87 GHz. 



Cummings model can describe the coupling of the qubit 
to the single mode cavity |5l|, [54j . 

In order to determine the conditions to obtain the 
model Hamiltonians we have used in the preceding sec- 
tions, thus we shall focus here the relatively less explored 
N-V center coupling to the cavity, in particular we exam- 
ine if such a coupling can be described by the Delta-type 
transition model. 

The diagram shown in Fig. [16] describe the relevant 
energy levels of a negatively charged N-V center [55|, [56j . 
The ground state manifold is 3 A, with zero field splitting 
D ~ 2.87 GHz, the excited state is 3 E separated from the 
ground level by the zero phonon line at 637 nm, and there 
is a metastable state x A with non-radiative transitions. 
Further splittings of 3 E level are in the order of 2.6 GHz 
and 2.3 GHz for 3 E X and 3 E yi respectively. 

Electric dipole (El) and magnetic dipole (Ml) selec- 
tion rules forbid an immediate realization of Raman type 
(A) transitions in the given energy level diagram, let 
alone a A transition. On the other hand, optical Ra- 
man induced spin-flipping processes J57H59J have been 
realized in addition to spin conserving (cyclic) transitions 
|60j in recent experiments. To comprehend how this hap- 
pens, detailed recent studies of the excited level reveal 
that among the potential spin-spin and spin-orbit inter- 
actions that can influence the excited level structure, it is 
the local strain induced level mixing that allows for the 
essential non-spin preserving transitions for an optical 
A-scheme. Using external magnetic or electric fields or 
specifically engineered N-V centers thus one can produce 
the level combinations for such transitions. To make our 
discussion concrete, we generalize the A-transition gen- 
eration protocol using external electric fields [61( to the 
case of A-transition. Typical electric field strengths sev- 
eral ~ MV/m are experimentally used |59j . 

Under strong electric field the level structure of the 
negatively charged N-V center is shown in the Fig. [T71 
Here, both cyclic (spin preserving) and spin-flipping tran- 
sitions can simultaneously happen. 
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FIG. 17. (Color online) Energy level diagram of negatively 
charged N-V center. ZPL is the zero phonon line at 638 nm 
(1.945 eV) and ZFS is the zero field splitting at 2.87 GHz. 



Optical fields are right circular polarized. El and ESR 
selection rules allows us to obtain cyclic transitions in 
this single N-V molecule in a hexagon scheme, where 
all nodes can be accessed starting from any other node, 
as shown in Fig. [T71 (b). The upper level A is an ex- 
cited level generated by the mixed states by the elec- 
tric field. The intermediate levels are \m s = — l,raj =t) 
and \m s = +l,raj =>|r), denoted by |— ) and |+), re- 
spectively. The lowest levels are \m s = 0,raj =t) and 
\m s = 0,772/ =>|r), denoted by |*f) and |^), respectively. 
Coupling of the cavity field to the ground state 3 A mani- 
fold is described by an anomalous Zeeman interaction of 
the form 



* = £*? 



Da 



fJ-B 



E3 



00* 



(31) 



where /i# is the electronic Bohr magneton and g is a diag- 
onal tensor with elements (g±, g±, g\\) with g± = 2.0024 
and g\\ = 2.0028 in a frame where the quantization 
axis (z) is chosen to be the symmetry axes of the N- 
V molecule. In the C$ v symmetry of diamond, one can 
find four of such axes corresponding to four crystallo- 
graphic classes NV — f,f = 1,2, 3, 4. We work a within 
a particular class, but omit its label (/) for now for 
the sake of notational simplicity. Expressions for dif- 
ferent classes can be obtained by using spherical tensor 
operations. Nuclear Zeeman term is too weak and ig- 
nored. Depending on the relative density of surrounding 
nuclei ( 1 4N, 1 5N, 1 3C), quadrupole and hyperfine inter- 
actions can impose significant decoherence on the cen- 
tral electronic spin [26|. Their contribution can be at- 
tributed to our effective loss terms in the master equa- 
tion treatment we used. The first term is the zero-field 
Hamiltonian, with D being the zero-field splitting ten- 
sor, which is diagonal in the molecular frame with the 
elements (D±,D±,D\\). It simplifies to Hzfs = Ba 2 z , 
up to a constant, with D = 3D\\/2 ~ 2.87 GHz. The 



index i = 1,...,A// is the index of N-V defects at loca- 
tion fi = (xi,yi,Zi) in the lab (cavity) frame. We shall 
consider a single crystal case first. 

Writing the cavity field as B(fi) = Bi(b + &t), with 
Bi = B_\_(xi,yi) cos (jrZi/X), the Hamiltonian becomes in 
the ground state basis 

H = Dj2Xa + 9B(b + tf)(X i+ + X t _ + /i.e.), (32) 
i=± 

where we neglected the small Zeeman shift, which is in 
the order of kHz for a cavity field in the order of nT, 
next to the D ~ GHz. Generalized Hubbard operators 
are introduced to be 



Xf 



X+- 



Here, B i± = B ix ± B ix and 



V^ B i+ R (i) 


(33) 


sr Bi ~ *?« 


(34) 


E*?- 


(35) 



s± = ,/£(*& + fl&). 



(36) 



In the large Nf and weak excitation limit, for a 
quasi-uniform magnetic field distribution, we can em- 
ploy the group contraction method as before, to get a 
bosonic representation for the Hubbard operators such 
that X±+ ~ B+ and X^_ ~ B_ . Consistently, they form 
the Schwinger representation of the remaining isospin 

subgroup so that we take X^ ~ B+B- which allows 

for replacing X ++ + X with B ] + B + + B ] _B_ in the 

zero field splitting term of the Hamiltonian. 

The interaction of the excited state with the optical 
fields can be described in terms of the collective operators 
as well. To make this consistently with the microwave 
transitions we first make the usual bosonization of ground 
level to excited level | A) transition via X^a ~ X^a ~ A. 
Here we implicitly assume the populations in degenerate 
ground levels are large and A-mode is weakly excited. 
Then consistency with the isospin groups is satisfied by 
taking their Schwinger representations as Xa± = A^B±. 
To achieve this limit adequately, we demand that the op- 
tical fields are to be made sufficiently homogeneous and 
isotropic over the diamond substances. This may be ac- 
complished by suitable confocal microscopy. The total 
Hamiltonian in this bosonic limit describes again a cyclic 
transition system which is reduced from a hexagon to a 
diamond scheme, as shown in Fig. [TT] (c). The degen- 
eracy of intermediate levels allows us to e mp l oy a sim- 
ple multilevel Morris- Shore transformation |62l-l64j in the 
form B = (B- + £+)/>/2 and B' = (B- - £+)/\/2. Ge- 
ometrically, this folds the diamond scheme into a delta- 
type cyclic transition as shown in Fig. [T71 (d) where B' 
is uncoupled from the system. 
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Even though we have described how to get a cyclic 
N-V ensemble, our analysis so far is limited to a single 
crystallographic class. To complete our discussion let us 
now address the case of all N-V classes. For that aim 
we make use of slightly more complicated Morris- Shore 
transformations . 

The Hamiltonian including all the N-V classes is 

4 4 

H = w A ^2 A \ A f + UB ^2 B \ B + ^a^a + u b b ] b 



f=i 



/=i 



UJQ 



-^a z + g(tf* + h.c.) + ft J2( A \ B f + h ' c ') 



/=i 



4 



(g B fbBl + #A/aA f f + ft.c.) 



We find that under the conditions 

#£2 #A2 #£3 #A3 #£4 #A4 



PB1 0A1 



<Ml PB1 #A1 



(37) 



(38) 



only two quasi-spin wave modes can be coupled to the 
electromagnetic modes. These so-called "bright" modes 
are given by collective bosonic modes 






where 



9a 



\ 



/=1 



9b 



\ 



5>B/- 



(39) 



(40) 



/=i 



The uncoupled, "dark" modes are listed in the appendix 
B. After that we arrive our starting Hamiltonians, Eqs. |4] 
-[5] in the Sec. HU A simple way to realize the dark quasi- 
spin wave modes is to make the interaction coefficients 
9a/ Bf equal for each N-V class. This is the case in the ex- 
periment of strong coupling of the N-V center and super- 
conducting cavity. By choosing a particular placement of 
the diamond crystal over the cavity substrate for which 
the cavity magnetic field makes effectively equal angles 
for the four quantization axes of the N-V classes, the 
Zeeman interaction coefficients become equal [26]. The 
collective enhancement of the interactions would depend 
on Nf, which we may assume Nf = N/4. Finally, the 
case of two diamond crystals can be immediately gener- 
alized from these results, where only the magnetic field 
would change a sign at the location of the other crystal 
half wavelength away. 

The typical parameters in these models are, according 
to recent experiments as discussed in Ref. [65j, as follows. 
Qubit decoherence times are in the order of few /is, in par- 
ticular for transmon is about ~ 4/iS [51] . For T\ ~ 7/xs 
and T2 ~ 500 ns, relaxation and dephasing rates are re- 
spectively 7i/2tt - 0.02 Mhz and 70/2tt - 0.31 MHz 
and the resonance frequency of the qubit is tunable in 
the ranges ujq ~ 5 — 15 GHz. Stripline cavity loss rate 



is k/2i\ = uj b /2i\Q ~ 5 kHz at uj b /2i\ = 5 GHz in a 
cavity with Q factor Q = 10 6 . Cavity frequency is ad- 
justable in the ranges usually about uj b /2ir = 5 — 10 GHz. 
Coupling strength between the cavity and the qubit can 
be adjusted in the ranges g ~ 5 — 200 MHz. In our 
numerical simulations we choose our parameters ranges 
compatible within these practical values. In particular 
we fix g a /Q = 0.1 so that J = — 0.1g&, E a = O.IQa- 
For g b = (2tt)11 MHz of Ref. (26J| , we take J/2tt = 1 
MHz. We assume Qa and detunings are adjusted to 
make E a = E b . E b = \/N cav K b with N cav = 0.01 is 
cavity photon number for the empty microwave cavity 
on probe resonance under weak drive conditions. We 
use 7 = 7i/2tt = 0.02 MHz and 7/2tt = 0.3 Mhz. We 
consider ranges of g/2ir up to 10 MHz and 5/2tt within 
(—5,5) MHz. To take into account the conditions of ex- 
periment Ref. [26| we use n b /2i\ = 0.4 MHz, while the 
effective loss of optical modes are chosen to be the same 
Ka — K>b-> which corresponds to ja/2tt ~ 40 MHz for the 
effective decay rate of the excited N-V center, which is 
reasonable with the experiments [66] . 

The quantum master equation we have used can be re- 
duced from a more general Bloch-Redfield master equa- 
tion [63, |68| under low temperature conditions [69] . 



V. CONCLUSION 

Summarizing, we have investigated quantum statisti- 
cal properties of optical radiation from a A-type cyclic 
atomic ensembles placed in a microwave cavity and cou- 
pled to a two- level atom. In particular, quantum co- 
herence and quantum correlations of the emitted optical 
photons are examined. Effective models describing the 
system analogous to coupled driven, dissipative linear 
and nonlinear resonators are obtained under Frochlich 
transformations . 

It is shown that analogous to the adiabatic transfer 
protocol between stationary qubits [28] , the quantum in- 
formation can be transferred directly between the sta- 
tionary qubit and the optical flying qubit. The proposed 
system completes the quantum memory and quantum 
hardware state transfer protocols by interfacing them 
with quantum communication channels. In addition, it 
is found that quantum coherence characteristics can be 
transferred between the microwave and optical fields. 

Furthermore, it is shown that distant entanglement be- 
tween two optical modes can be realized in the case of 
two-cyclic atomic ensembles. In addition conditions of 
three- mode entanglement are also revealed. Controllable 
particle and mode entanglement for the optical modes or 
between the optical and the microwave photons is shown 
to be realizable on demand. 

Finally, diamond crystals with N-V centers in super- 
conducting coplanar waveguide resonators coupled to 
transmon qubits are proposed to physically realize the 
model system. The parameters used in numerical sim- 
ulations are justified for the proposed system. Under 
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strong electric field, it is described that a single N-V cen- 
ter allows for a particular hexagon- type cyclic transitions. 
Using group contraction method for the ensemble case, 
the transition type reduces to a diamond- type. For a gen- 
eral treatment of the N-V center taking into account all 
four crystallographic classes, and using generalized multi- 
level Morris- Shore transformations a further reduction to 
A-type cyclic scheme is obtained. 

We hope our work can inspire and contribute ongo- 
ing efforts for interfacing of optical networks to quantum 
memories and information processing units. 
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Appendix A: Coefficients of Effective Hamiltonian 



Here we list the coefficients in Eq. [T2l explicitly. 



In the ft ^> A limit, these coefficients can be approxi- 
mated to be 






J: 



2 ft 2 ' - " ' 4 

A 2 



1 



(At + A/2) 2 + 8(A a - A/2) 2 
2ft 2 



(A10) 

(All) 

(A12) 



0^^ + |(A 6 + ^) + |(A a -^A13) 



Q ~ 2ft 2 



G±«± 

E a =- 



99b 

y/2tt 

&Agb 
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E' b = E b 
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where A, \i = ±. 
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Appendix B: Dark quasi-spin wave modes 

Here we list the dark quasi-spin wave modes that are 
uncoupled from the N-V ensemble Hamiltonian 



B' 



B" 



gB34 

gBgBi2 
gBi2 



gBgB34 



■(gBiBi J cgB2B 2 ) 

(#£3^3 + gBAB^) 



gB2 
gB\2 



B x 



B"> = 1^B 3 



gB34 



gm 

gBi2 
gB3 

gB34 



B 2 



B±, 



where we define 



g%i 



gB2-> 



gB34 



gm 



g%4- 



The same definitions and dark modes for the 
can be simply obtained by replacing B — >• A 
equations. 
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